Slow dynamics at the smeared phase transition of randomly layered magnets 
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We investigate a model for randomly layered magnets, viz. a three-dimensional Ising model with 
planar defects. The magnetic phase transition in this system is smeared because static long-range 
order can develop on isolated rare spatial regions. Here, we report large-scale kinetic Monte Carlo 
simulations of the dynamical behavior close to the smeared phase transition which we characterize 
by the spin (time) autocorrelation function. In the paramagnetic phase, its behavior is dominated 
by Griffiths effects similar to those in magnets with point defects. In the tail region of the smeared 
transition the dynamics is even slower: the autocorrelation function decays like a stretched expo- 
nential at intermediate times before approaching the exponentially small asymptotic value following 
a power law at late times. Our Monte-Carlo results are in good agreement with recent theoretical 
predictions based on optimal fluctuation theory. 

PACS numbers: 64.60.Ht, 05.50.+q, 75.10.Nr, 75.40.Gb 



I. INTRODUCTION 

In recent years, there has been a resurgent interest 
in the influence of defects, impurities or other types 
of quenched disorder on the properties of phase transi- 
tions and critical phenomena. This renewed attention is 
largely due to the discoveries of novel disorder effects that 
go beyond the framework of perturbation theory and the 
perturbative renormalization group. 

A particularly interesting class of non-perturbative 
phenomena are the so-called Griffiths or rare region ef- 
fects that are produced by rare strong spatial disorder 
fluctuations. They can be easily understood on the ex- 
ample of a diluted magnet: Due to the dilution, the crit- 
ical temperature T c of the disordered system is lower 
than its clean value, T . In the temperature interval 
T c < T < T c °, the diluted bulk system is in the disor- 
dered phase. However, in an infinite size sample, there is 
an exponentially small, but nonzero probability for find- 
ing an arbitrary large spatial region devoid of impurities. 
Such a rare region, a 'Griffiths island', displays local or- 
der in the temperature interval T c < T < T®. Due to its 
size, such an island will have very slow dynamics because 
flipping it requires changing the order parameter over a 
large volume. More than 30 years ago, Griffiths 1 showed 
that the presence of these locally ordered islands pro- 
duces a singularity in the free energy in the whole region 
T c < T < T c °, which is now known as the Griffiths region 
or the Griffiths phase. 2 In generic classical systems with 
uncorrelated or short-range correlated disorder, thermo- 
dynamic Griffiths effects are very weak because the sin- 
gularity in the free energy is only an essential one. 3-5 
To the best of our knowledge, classical thermodynamic 
Griffiths singularities have therefore not been verified in 
experiment (see also Ref. 6). 

In contrast to the thermodynamics, the long-time dy- 
namics is dominated by the rare regions. Inside the Grif- 
fiths phase, the spin autocorrelation function C(t) decays 
very slowly with time t, as lnC(t) ~ — (lni) d /( d_1 ) for 
Ising systems 2,7-10 and as In C(t) ~ — t 1 / 2 for Heisen- 



berg systems. 9,11 These results were recently confirmed 
by more rigorous calculation for the equilibrium 12,13 and 
dynamic 14 ' 15 properties of disordered Ising systems. 

The effects of impurities and defects are greatly en- 
hanced by long-range spatial disorder correlations. In 
particular, if the disorder is perfectly correlated in some 
spatial directions, the rare regions are extended objects 
which are infinite in the correlated dimensions. This 
makes their dynamics even slower and so increases their 
effects. In Ising models with linear defects (disorder per- 
fectly correlated in one dimension), the thermodynamic 
Griffiths singularities are of power-law type, with the av- 
erage susceptibility actually diverging in a finite tem- 
perature region. The critical point itself is an exotic 
infinite-randomness critical point and displays activated 
rather than power-law scaling. This was first found in 
the McCoy- Wu model, 16 ' 17 a disordered 2D Ising model 
in which the disorder is perfectly correlated in one dimen- 
sion and uncorrelated in the other. Later it was studied 
in great detail in the context of the quantum phase tran- 
sition of the random transverse-field Ising model where 
the imaginary time dimension plays the role of the "cor- 
related" direction. 18 ' 19 

Recently, it has been shown that even stronger effects 
than these power-law Griffiths singularities can occur in 
Ising models with planar defects. 20 ' 21 Because the disor- 
der is perfectly correlated in two directions, the effective 
dimensionality of the rare regions is two. Therefore, an 
isolated rare region can undergo the phase transition in- 
dependently from the bulk system. This leads to a de- 
struction of the global sharp phase transition by smear- 
ing. Similar smearing effects have also been found in 
itinerant quantum magnets 22 and at a non-equilibrium 
transition in the presence of linear defects 23 . A recent 
review of these and other rare region effects can be found 
in Ref. 24. 

In this paper, we study the dynamics of an Ising model 
with planar defects in the vicinity of this smeared phase 
transition by large-scale kinetic Monte-Carlo simulations. 
The paper is organized as follows. In Section II we in- 
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troduce the model and briefly summarize the results of 
the optimal fluctuation theory 20,24 ' 25 for smeared phase 
transitions to the extent necessary for the analysis of our 
data. In section III, we explain the simulation technique, 
we present our results for the spin autocorrelation func- 
tion, and we compare them to the theoretical predictions. 
Conclusions are presented in section IV. 



II. THEORY 

A. Three-dimensional Ising model with planar 
defects 

We consider a three-dimensional classical Ising fer- 
romagnet with planar defects, the same model whose 
thermodynamics was studied in Ref. 21. Ising variables 
Sij^k — ±1 reside on the sites of a cubic lattice. They 
interact via nearest-neighbor interactions. In the clean 
system all interactions are identical and have the value J. 
The defects are modeled via 'weak' bonds randomly dis- 
tributed in one dimension (uncorrelated direction). The 
bonds in the remaining two dimensions (correlated direc- 
tions) remain equal to J. The system effectively consists 
of strongly-coupled slabs (layers) of varying thickness, 
separated by parallel layers of weak bonds. The Hamil- 
tonian of the system is given by: 

H = — JiSij^Si+ij.k 

i=l,...,L ± 
J,k = l,...,L c 

y ' J(Si t j^Si t j-\-\ t k + <S'i, i j,fc<S'i, i j,fc+l))(l) 



i=l,...,L ± 
j,k = l,...,L c 

where L±(Lc) is the length in the uncorrelated (corre- 
lated) direction, i, j and k are integers counting the sites 
of the cubic lattice, J is the interaction energy in the cor- 
related directions and Ji is the random interaction in the 
uncorrelated direction. The Ji are drawn from a binary 
distribution: 



Ji = 



cJ with probability p 
J with probability 1 — p 



(2) 



characterized by the concentration p with < p < 1 and 
the relative strength c of the weak bonds (0 < c < 1). 

The Hamiltonian (1) does not contain any internal dy- 
namics. Instead, the dynamics must be added ad-hoc 
by specifying an appropriate dynamical algorithm. We 
consider a purely relaxational local dynamics without 
any conservation laws, i.e., model A in the Hohcnbcrg- 
Halperin classification. 26 Microscopically, it can be real- 
ized, e.g., by the Glauber or Metropolis algorithms. 27 ' 28 



B. Smeared phase transition and optimal 
fluctuation theory 

In the absence of impurities (p = 0), the system (1) 
undergoes a ferromagnetic phase transition at the clean 



critical temperature T c ° = 4.511J, with the order param- 
eter being the total magnetization: 



m 



(3) 



Here, V = L±L C is the volume of the system, and (•) 
denotes the thermodynamic average. 

In the presence of disorder (p > 0), a crucial role is 
played by rare strong disorder fluctuations: In analogy 
to the Griffiths phenomena 1 discussed in the Introduc- 
tion, there is a small but finite probability for finding a 
large spatial region containing only strong bonds in the 
uncorrelated direction. Such a rare region can be locally 
in the ferromagnetic state while the bulk system is still 
in the disordered (paramagnetic) phase. The ferromag- 
netic order on the largest rare regions starts to emerge 
right below the clean critical temperature T°. Since the 
defects in the system are planar, these rare regions are 
infinite in the two correlated dimensions but finite in the 
uncorrelated direction. Each rare region is thus equiva- 
lent to a two-dimensional Ising system that can undergo 
a real phase transition independently of the rest of the 
system. The resulting effect is much stronger than con- 
ventional Griffiths effects: the global phase transition is 
destroyed by smearing, and the order parameter develops 
very inhomogeneously in space with different parts of the 
system (different rare regions) ordering independently at 
different temperatures. 20,21,24 

The leading thermodynamic behavior in the tail of 
the smeared transition can be determined using opti- 
mal fluctuation theory. 20,24 The approach is similar to 
that of Lifshitz 29 and others for the description of the 
tails in the electronic density of states of disordered sys- 
tems. The theory can be easily developed for a general 
d-dimensional system with dc correlated dimensions and 
d± = d — dc uncorrelated dimensions. In the case of the 
Hamiltonian (1), dc — 2 and d± = 1. 

In the tail region of the smeared transition, the sys- 
tem consists of a few isolated ferromagnetic rare regions 
embedded in a paramagnetic bulk. We start from the 
probability w for finding a large region of linear size L RR 
containing only strong bonds. Up to pre-exponential fac- 
tors, it is given by 



w ~ (1 — p) 1 



cxp 



4Vn(i-p) 



(4) 



Such a rare region develops static ferromagnetic long- 
range order at some temperature T C (L RR ) below the 
clean critical temperature T°. The value of T C (L RR ) 
varies with the size of the rare region; the largest is- 
lands will develop long-rage order closest to the clean 
critical point. To determine T C (L RR ), we can use finite 
size scaling 30 for the clean system because a rare region 
is equivalent to a "slab" of a clean Ising model. This 
yields 



T°-T C (L RR ) = \r c (L RR )\ 



A J-V 



(5) 
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where <f> is the finite-size scaling shift exponent of the 
clean system and A is the amplitude for the crossover 
from d = dc + d± dimensions to a slab geometry infi- 
nite in dc (correlated) dimension but with finite length 
in the other (uncorrelated) directions. The reduced tem- 
perature r = T — Xj? measures the distance from the clean 
critical point. As long as the clean d-dimensional Ising 
model is below its upper critical dimension (d+ = 4), 
hyperscaling is valid and the finite-size shift exponent is 
related to the correlation length exponent v by (f> = 1/v 
which we assume from now on. Combining (4) and (5), 
we obtain the probability for finding an island which be- 
comes critical at some r c as: 



w(r c ) ~ exp (-B\r c \- d± ") (for r. 



0- 



(6) 



with the constant B = - ln(l - p)A d±l/ . The total (av- 
erage) magnetization m at some reduced temperature r 
is obtained by integrating over all rare regions which are 
ordered at r, i.e., those with r c > r. Since the functional 
dependence on r of the local magnetization on the island 
is of power-law type it does not enter the leading expo- 
nentials but only pre-exponential factors. To exponential 
accuracy, we therefore obtain 

m(r) ~cxv(-B\r\- d± ») (for r -> 0-). (7) 

Thus, the total magnetization develops an exponential 
tail towards the disordered phase which reaches all the 
way to the clean critical point. Analogous estimates show 
that the homogeneous magnetic susceptibility does not 
diverge anywhere in the tail region of the smeared tran- 
sition. At r = 0, the exponentially decreasing island 
density overcomes the power-law divergence of the sus- 
ceptibility of an individual island; and once a nonzero 
magnetization has developed, it cuts off any possible sin- 
gularity. However, there is an essential singularity in the 
free energy at the clean critical point produced by the 
vanishing density of ordered islands. 



C. Dynamics at the smeared transition 

After having briefly discussed the optimal fluctuation 
theory for the thermodynamics, we now consider the dy- 
namical behavior at the smeared phase transition of our 
disordered magnet (1). The interesting physics in the tail 
of the smeared transition is local with respect to the un- 
correlated dimensions because different rare regions are 
effectively decoupled from each other, and the spatial 
correlation length in these directions remains finite. An 
appropriate quantity to study the rare region dynamics 
is therefore the time autocorrelation function of the Ising 
spins, 



(8) 



The leading long-time behavior of C{t) in the tail of the 
smeared transition can be determined using optimal fluc- 
tuation arguments similar to that of section II B. 24,25 



According to finite-size scaling 30 , the behavior of the 
correlation time £ t of a single rare region of size L RR 
in the vicinity of the clean bulk critical point can be 
modelled by (for T < T c °, i.e., r < 0) 



it{r,L RR ) ~ L 



{zv — zD) I v 
RR 



AL 



RR 



(9) 



Here, v and z are the correlation length and dynamical 
exponents of a d-dimensional system, and v and z are the 
corresponding exponents of a dc-dimensional system. 

Let us first consider the time evolution of the auto- 
correlation function C(t) at the clean critical point T°, 
i.e., the boundary between the conventional paramag- 
netic phase and the smeared tail of the ordered phase. 
For r = 0, the correlation time (9) simplifies to £ t ~ L RR . 
The rare region contribution to C(t) is obtained by sim- 
ply summing over the exponential time dependencies of 
the individual islands with L^-dependent correlation 
time. Using (4) we obtain to exponential accuracy 



C(t) ~ J dL RR exp (4V*(1 ~P) - Dt/L RR ) 



(10) 



where D is a constant. This integral can be easily es- 
timated within the saddle-point method. The leading 
long-time decay of the autocorrelation function C(t) at 
the clean critical point is given by a stretched exponen- 
tial, 



lnC(i) ~ -[-ln(l -p)y-/( d ^+*) t d±/(d±+z) ( r = 0). 

(11) 

In the conventional paramagnetic phase, T > T c ° 
(r > 0), the correlation time does not diverge for any 
L RR . Instead, the correlation time of the large islands 
saturates at £ t {r,L RR ) - r- zv for L RR > {r/A')- u . The 
autocorrelation function C(t) can again be evaluated as 
an integral over all island contributions. We find that 
there are two regimes separated by a crossover time 



-{d ± +z)v 



(12) 



For intermediate times t < t x , the autocorrelation func- 
tion follows the stretched exponential (11). For times 
larger than the crossover time t x , we obtain a simple ex- 
ponential decay 

lnC(i)~-Vr (for t > t x , r > 0) (13) 

with the decay time r ~ r~ zu . Our results for T > T 
agree with the corresponding conventional dynamical 
Griffiths effects 9 for systems with point defects. This 
is not surprising, because above Tj?, there is no quali- 
tative difference between the Griffiths and the smearing 
scenarios: In both cases, all rare regions are locally still 
in the disordered phase. 

This changes in the tail of the smeared transition below 
the clean critical point Tj . For r < 0, we repeat the 
saddle point analysis with the full expression (9) for the 
correlation length. Again, for intermediate times t <t x , 
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the decay of the average density is given by the stretched 
exponential (11). For times larger than the crossover 
time t x the system realizes that some of the rare regions 
have developed static order and contribute to a non-zero 
long-time limit of the autocorrelation function C(t). The 
approach of C(t) to this limiting value is dominated by 
finite-size islands with Lrr ~ (— rjA)~ v because they 
have diverging correlation time. As a result, we obtain a 
power-law, 



C(t) - C(oo) ~ r 



(for t>t x , r < 0). (14) 



The value of ip cannot be found by our methods since it 
depends on the neglected pre-exponential factors. 

In Ref. 25, the predictions (11), (13), and (14) have 
been compared to results from a local dynamical mean- 
field theory. The purpose of the present paper is to verify 
these relations by large-scale kinetic Monte-Carlo simu- 
lations of a realistic model with short-range interactions, 
viz. the Hamiltonian (1). 




FIG. 1: Total magnetization m and susceptibility \ as func- 
tions of T for L± — 100, Lc = 200 and p — 0.2 averaged over 
200 disorder realizations (from Ref. 21). 



III. COMPUTER SIMULATIONS 
A. Method 

We now turn to the main part of the paper, kinetic 
Monte-Carlo simulations of our three-dimensional Ising 
model with nearest neighbor interactions and planar 
bond defects, as given in cq. (1). Since the smearing 
of the transition is a result of exponentially rare events, 
sufficiently large system sizes are required in order to ob- 
serve it. In the production runs, we have mostly used 
system sizes of 260 3 or 300 3 sites. In this context, let 
us briefly comment on finite-size effects. Because the 
phase transition in an Ising model with planar defects 
is smeared, conventional finite-size scaling 30 does not di- 
rectly apply. Nonetheless, there are two types of finite- 
size effects, whose consequences for the thermodynamics 
were discussed in detail in Refs. 20,21. A finite size in 
the uncorrelated direction limits the possible sizes of the 
rare regions and thus cuts off the probability distribution 
(6), while a finite size in the correlated directions rounds 
the local phase transitions of the rare regions. We have 
carried out test runs using 50 3 , 100 3 and 200 3 sites to 
ensure that our results are not influenced by these finite- 
size effects. 

We have chosen J = 1 and c = 0.1 in cq. (2), i.e., 
the strength of a 'weak' bond is 10% of the strength of 
a strong bond. The simulations have been performed for 
various disorder concentrations p = {0.2, 0.3, 0.4, 0.5, 0.6} 
over a temperature range of T = 4.25 to 5.00, covering 
both sides of the clean critical point at T c ° = 4.511. To 
suppress the statistical disorder fluctuations we have av- 
eraged the results over a large number of disorder config- 
urations ranging from 100 to 5000 depending on param- 
eters. 

For each disorder configuration, the simulation pro- 
ceeds as follows: After the initial setup, the system 



is equilibrated using the highly efficient Wolff cluster 
algorithm. 31 (Using a cluster algorithm is possible for 
our Hamiltonian because the disorder does not induce 
frustration.) The length of the equilibration period is 50 
Wolff sweeps. (A Wolff sweep is defined as a number of 
cluster flips such that the total number of flipped spins 
is equal to the number of sites, i.e., on the average each 
spin is flipped once per sweep.) The actual equilibration 
times have typically been of the order of 10-20 sweeps at 
maximum. Thus, an equilibration period of 50 sweeps is 
more than sufficient. 

After the system is equilibrated, the spin configuration 
is stored and the simulation time is set to t = 0. The ac- 
tual time evolution of the system according to model A 
dynamics is performed using the Metropolis algorithm. 27 
After each sweep (a Metropolis sweep consists in one at- 
tempted spin flip per site), the spin configuration is deter- 
mined and the spin autocorrelation function is calculated 
from (8). The length of this measurement period is up 
to 5000 Monte-Carlo (Metropolis) sweeps. 



B. Results 

In this section, we present the results of our kinetic 
Monte-Carlo simulations and compare them to the the- 
oretical predictions of section II C. For later reference, 
we first give a brief overview of the thermodynamics, as 
determined in Ref. 21. Fig. 1 shows total magnetization 
and susceptibility as functions of temperature averaged 
over 200 samples of size L± = 100 and Lc = 200 with an 
impurity concentration p = 0.2. For temperatures above 
T = 4. 2... 4. 3, the magnetization shows a pronounced 
tail, i.e., it vanishes very gradually when approaching 
the clean critical temperature T = 4.511. It should be 
emphasized that this tail is not due to finite-size effects. 
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FIG. 2: Spin autocorrelation function C as a function of time 
t for L± = Lc = 260, p — 0.3. The number of disorder 
realizations ranges from 100 at the lowest temperatures to 
5000 at the highest temperatures. The power-law time axis 
corresponds to the stretched exponential predicted in (11). 

Indeed, a detailed analysis 21 has shown that the magne- 
tization nicely follows the prediction (7) of the smeared 
transition scenario over more then a magnitude in m. 
Analogously, an analysis of the susceptibility has shown 
that it does not diverge in the thermodynamic limit, in- 
stead, it displays a rounded peak. 

We now turn to the spin autocorrelation function C'(t). 
Figure 2 gives an overview over its behavior for an im- 
purity concentration p = 0.3 and times up to t = 5000 
over the entire temperature range studied. In this figure 
the data are plotted in the form InC vs. t 9 - 329 such that 
the stretched exponential (11) predicted for the time de- 
pendence at the clean critical temperature T 9 = 4.511 
gives a straight line. Here we have used d± = 1 and a 
value of z = 2.04 for the dynamical exponent of the clean 
three-dimensional Ising model with model A dynamics 
(see, e.g., Refs. 32-34). The figure shows that the au- 
tocorrelation function at T 9 indeed follows the stretched 
exponential over more than two orders of magnitude in 
C. At temperatures T ^ T 9 , the autocorrelation func- 
tion initially follows the same stretched exponential but 
eventually crosses over to a different asymptotic form: 
For T > T 9 , it decays faster than the stretched expo- 
nential, while it seems to saturate at a nonzero value for 
T < T 9 . In the following, we investigate the different 
temperature regions in more detail. 

Let us begin by analyzing the autocorrelation func- 
tion C(t) right at the clean critical temperature T®. 
Figure 3 shows C(t) for several impurity concentrations 
p = 0.2... 0.6, again plotted as InC vs. i 329 . After 
the initial transients, all curves are straight lines, indi- 
cating that C(t) follows the theoretical prediction (11) 
over at least two orders of magnitude in C for all concen- 
trations. From the curves in Fig. 3, one can determine 
the decay constants (i.e., the slopes) as a function of im- 
purity concentration p. According to (11), the decay con- 




r 0.329 



FIG. 3: Spin autocorrelation function C(t) at the clean crit- 
ical temperature T<? — 4.511 for L± = Lc = 300 and several 
impurity concentrations p. The data are averages over 100 
disorder realizations (1000 realizations for p — 0.3). The sta- 
tistical errors are not bigger than a symbol size for C > 10~ 3 . 
For p = 0.3, they reach about twice the symbol size at the 
right end of the curve. 




[-In(1-P)f 671 

FIG. 4: Decays constants (slopes) of the stretched exponen- 
tials of Fig. 3 as a function of impurity concentration p, plot- 
ted according to the prediction (11). 

stant should be proportional to [— ln(l — p)] z /( d ^+ z ) = 
[— ln(l — p)] 671 . Figure 4 shows the corresponding plot 
for our data. It shows that the behavior of the decay 
constants is in good agreement with the theoretical pre- 
dictions. 

After the behavior at the clean critical temperature T 9 , 
we now consider the conventional paramagnetic phase 
T > T 9 . In Figure 5, we plot the autocorrelation func- 
tion C(t) for p = 0.3, L± = Lq = 260, and temperatures 
T > T 9 as InC vs t such that a simple exponential de- 
cay gives a straight line. For all shown temperatures, 
the long-time behavior of the autocorrelation function is 
indeed an exponential decay in agreement with the the- 
oretical prediction (13). To determine the decay time r 
we fit the long-time behavior of C(t) to (13). The in- 
set of Fig. 5 shows r _1 as a function of (T — T 9 ) 2I/ with 
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FIG. 5: Log-linear plot of the spin autocorrelation function for 
temperatures above T° (from the top: T = 4.6,4.7,4.8,5.0). 
The other parameters are as in Fig. 2. The solid lines are fits 
of the long-time behavior to (13). The statistical errors are 
not bigger than a symbol size for C > 10 -4 . Inset: Inverse 
decay time r _1 as a function of (T — T°) zv . 




t 



FIG. 6: Double- logarithmic plot of the dynamical part C(t) — 
C(oo) of the autocorrelation function for temperatures below 
T°. The other parameters are as in Fig. 2. The straight 
lines are fits of the long-time behavior to the power law (14) 
giving exponents of if) = 0.80 (T = 4.35), 0.90 (T = 4.40) and 
1.23 (T = 4.45). 



the three-dimensional clean correlation length exponent 
given by v = 0.628 (see, e.g., Ref. 35) and z = 2.04 as 
before. As predicted by the optimal fluctuation theory 
in Sec. II C, the inverse decay time depends linearly on 
(T — T®) zv . (The remaining small deviations can proba- 
bly be attributed to the pre-exponential factors neglected 
in the optimal fluctuation theory.) We have performed 
analogous simulations for an impurity concentration of 
p = 0.5 and found the same qualitative behavior. Of 
course, non-universal prefactors have different values. 

Finally, we turn to the properties of the spin autocorre- 
lation function C (t) below the clean critical temperature 
T c °, i.e., in the tail of the smeared phase transition. For 
these temperatures, the total equilibrium magnetization 
is nonzero because some of the rare regions have already 
developed a static magnetization (see also Fig. I). 20 ' 21 
Consequently, the autocorrelation function does not de- 
cay to zero in the long-time limit but rather approaches 
the limiting value m 2 as can be seen in the overview fig- 
ure 2. The dynamical correlations of the magnetization 
fluctuations (Sij t k — m) are represented by the devia- 
tion C(t) — C(oo) of the autocorrelation function from 
its long-time limit. Figure 6 shows a double-logarithmic 
plot of C{t) - C*(oo) for p = 0.3, L± = L c = 260, and 
temperatures T = 4.35, 4.40, and 4.45. For all three tem- 
peratures, the long-time behavior of C(t) — C(oo) follows 
a power law. Fits to (14) give exponents of if) = 0.80 
for T = 4.35, 0.90 for T = 4.40, and 1.23 for T = 4.45. 
Thus, as predicted, the exponent of the power-law (14) is 
nonuniversal. (We note that the apparent crossing of the 
curves in Fig. 6 is not related to the long-time dynamics 
but simply the result of the subtraction of C(oc). At 
lower temperatures, C(oo) is larger, and since C(t) can- 
not exceed 1, the difference C(t) — C(oo) at early times 
must decrease with decreasing T. Since in the long-time 



limit, the decay is faster for larger T, the curves must 
cross.) 

We have carried out analogous simulations for the im- 
purity concentration p = 0.5, and we have obtained 
equivalent results. 

IV. CONCLUSIONS 

To summarize, we have studied the dynamic behav- 
ior of randomly layered Ising magnets by performing 
large-scale kinetic Monte-Carlo simulations of a three- 
dimensional Ising model with planar defects. In this sys- 
tem, the magnetic phase transition is smeared because 
rare strongly coupled spatial regions (i.e., thick strongly- 
coupled layers) independently undergo the phase tran- 
sition. We have found that the dynamics in the rare- 
region dominated tail of the smeared transition is very 
slow. The spin autocorrelation function approaches its 
stationary value following a power-law in time. At the 
clean critical temperature T~ (which marks the bound- 
ary between the conventional paramagnetic phase and 
the tail of the smeared transition), the autocorrelation 
function decays like a stretched exponential in time. In 
the paramagnetic phase above Tj?, the decay is of simple 
exponential type. However, the decay time r of this ex- 
ponential diverges in the limit T — > T®+. Our numerical 
results (both the functional forms and the exponent val- 
ues of the various time and disorder dependencies) are in 
very good agreement with a recent optimal fluctuation 
theory. 25 

All our explicit results are for strong impurities (c = 
0.1). We have chosen this value because strong impu- 
rities allow us to easily observe the smeared transition 
in a finite-size simulation. If the impurities are weak (c 
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close to 1), the smeared transition is too close to the 
clean critical point. Strong bulk fluctuations thus com- 
pete with the smearing, and the latter can only be ob- 
served in larger systems and/or at longer times. The 
behavior of our model for c = (dilution) is special, be- 
cause for this value, the system is decomposed into an 
ensemble of noninteracting slabs. 

We emphasize that we have considered a purely re- 
laxational (local) dynamics corresponding to model A in 
the Hohenberg-Halpcrin classification. 26 Other dynamic 
algorithms require separate investigations. For instance, 
model B dynamics globally conserves the order parame- 
ter. In this context, an interesting question is: How does 
the interplay of the local thermodynamics of the rare re- 
gions and the global conservation law modify the dynamic 
behavior in the tail of the smeared transition. 

Let us briefly compare the dynamics at a smeared 
phase transition studied here to the behavior in a con- 
ventional Griffiths phase (as produced by uncorrelated 
or short-range correlated disorder such as point defects). 
Both in a Griffiths phase and at a smeared transition, 
rare regions dominate the long-time dynamics. In a con- 
ventional Griffiths phase, the finite-size rare regions re- 
main fluctuating for all temperatures above the dirty 
critical point, i.e., their island correlation times remain 
finite. As a result, the autocorrelation function decays 

like lnC(t) (lnt) d ^ d -^ for Ising systems 2 - 7 " 10 and 

lnC(t) ~ — t 1 / 2 for Heisenberg systems 9 ' 11 . In contrast, 
in the tail of a smeared transition, the effects of the rare 
regions are even stronger because individual islands can 
undergo the phase transition independently, connected 
with a divergent island correlation time. Summing over 



all islands then leads to an even slower power-law decay 
of the spin autocorrelation function. 

Let us also comment on experiments. A direct realiza- 
tion of the scenario discussed here could be achieved by 
growing alternating layers of two ferromagnetic materials 
with different critical temperatures T c . To introduce dis- 
order, the thickness of the layers should be random. In 
the tail of the smeared transition, i.e., for temperatures 
close to but below the higher of the two T c , the dynamic 
magnetic response of this system will be dominated by 
the rare region contributions. Specifically, the dynamic 
susceptibility x( w ) wm be dominated by its local part 
which is essentially the Fourier transform of the autocor- 
relation function. Thus measuring the dynamic response 
in such a randomly layered magnet should provide an 
experimental verification of our results. 

Finally, while the results here have been derived for 
an Ising model with planar defects, we expect analogous 
results for other disorder-smeared phase transitions. In- 
deed, in the tail of the smeared non-equilibrium phase 
transition of a contact process with extended impurities, 
a power-law decay of the density was recently found. 23,36 . 
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